Source of Acquisition 
NASA Goddard Space Flight Center 


Precision Fe K a and Fe K (3 Line Spectroscopy of the 
Seyfert 1.9 Galaxy NGC 2992 with Suzaku 

Tahir Yaqoob , 1 ’ 2 Kendrah D. Murphy , 1 Richard E. Griffiths , 3 Yoshito Haba , 4 
Hajime Inoue , 5 Takeshi Itoh , 4 Richard Kelley , 2 Motohide Kokubun , 6 
Alex Markowitz , 2 ’ 7 Richard Mushotzky , 2 Takashi Okajima , 1 ’ 2 Andrew Ptak , 1 ’ 2 
James Reeves , 1,2 Peter J. Serlemitsos , 2 Tadayuki Takahashi , 5 and Yuichi Terashima 5 
1 Department of Physics and Astronomy, Johns Hopkins University, 

3400 N. Charles St., Baltimore, MD 21218 
2 Exploration of the Universe Division, NASA Goddard Space Flight Center, 

Greenbelt Rd., Greenbelt, MD 20771 

3 Department of Physics, Carnegie Mellon University, 5000 Forbes Avenue, 

Pittsburgh, PA 15213 

4 Department of Astrophysics, School of Science, Nagoya University, Chikusa-ku, 

Nagoya 464-01, Japan 

5 Institute of Space and Astronautical Science, 

Japan Aerospace Exploration Agency, 3-1-1 Yoshino-dai, 

Sagamihara, Kanagawa 229-8510, Japan 
6 Department of Physics, University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, 

Tokyo, Japan 

7 NASA Postdoctoral Research Associate 

(Received DRAFT VERSION 30 July 2006; accepted TBD 2006) 

Abstract 

We present detailed time-averaged X-ray spectroscopy in the 0.5-10 keV band 
of the Seyfert 1.9 galaxy NGC 2992 with the Suzaku X-ray Imaging Spectrometers 
(XIS). The source had a factor ~ 3 higher 2-10 keV flux (~ 1.2 x 10 -11 erg cm -2 s _1 ) 
than the historical minimum and a factor ~ 7 less than the historical maximum. 

The XIS spectrum of NGC 2992 can be described by several components. There is a 
primary continuum, modeled as a power-law with a photon index of P = 1.57lo!o3 that 
is obscured by a Compton-thin absorber with a column density of 8.0loJ x 10 21 cm -2 . 
There is another, weaker, unabsorbed power-law component (modeled with the same 
slope as the primary), that is likely to be due to the primary continuum being electron- 
scattered into our line-of-sight by a region extended on a scale of hundreds of parsecs. 

We measure the Thomson depth of the scattering zone to be r es = 0.072 ±0.021. An 
optically-thin thermal continuum emission component, which probably originates in 
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the same extended region, is included in the model and yields a temperature and 
luminosity of kT = 0.656lo!o6i keV and ~ 1.2 ±0.4 x 10 40 erg s _1 respectively. We 
detect an Fe K emission complex which we model with broad and narrow lines and 
we show that the intensities of the two components are decoupled at a confidence 
level > 3a. The broad Fe Ka line has an equivalent width of 1 eV and could 
originate in an accretion disk (with inclination angle greater than ~ 30°) around 
the putative central black hole. The narrow Fe Kct line has an equivalent width of 
163lfg eV an d is unresolved (FWHM < 4630 km s _1 ) and likely originates in distant 
matter. The absolute flux in the narrow line implies that the column density out of 
the line-of-sight could be much higher than measured in the line-of-sight, and that 
the mean (historically-averaged) continuum luminosity responsible for forming the 
line could be a factor of several higher than that measured from the data. We also 
detect the Fe K/3 line (corresponding to the narrow Fe Ka line) with a high signal- 
to-noise ratio and describe a new robust method to constrain the ionization state of 
Fe responsible for the Fe Ka and Fe K f3 lines that does not require any knowledge of 
possible gravitational and Doppler energy shifts affecting the line energies. For the 
distant line-emitting matter (e.g. the putative obscuring torus) we deduce that the 
predominant ionization state is lower than Fe VIII (at 99% confidence) , conservatively 
taking into account residual calibration uncertainties in the XIS energy scale and 
theoretical and experimental uncertainties in the Fe K fluorescent line energies. From 
the limits on a possible Compton-reflection continuum it is likely that the narrow 
Fe Ka and Fe K/3 lines originate in a Compton-thin structure. 

Key words: galaxies: active - galaxies: Seyfert - line:profiles - X-rays:galaxies 
- X-rays: individual (NGC 2992) 


1. Introduction 

The Fe K emission line profiles in active galactic nuclei (AGN) generally consist of a 
prominent core centered at ~ 6. 4 keV (typically with FWHM < 10,000 km s 1 or so) and 
in some sources an additional underlying broader component. The relatively narrow core at 
~ 6.4 keV is found to be ubiquitous in nearby, low to moderate luminosity AGN (e.g. Nandra 
et al. 1997a, b; Sulentic et al. 1998; Weaver, Gelbord, & Yaqoob 2001; Reeves 2003; Page 
et al. 2004; Yaqoob &; Padmanabhan 2004; Jimenez-Bailon et al. 2005; Zhou & Wang 2005; 
Jiang, Wang, & Wang 2006; Levenson et al. 2006). The broad Fe K line component potentially 
carries signatures of the inner accretion disk in the vicinity of the putative black hole by virtue 
of gravitational and Doppler energy shifts (e.g. Tanaka et al. 1995; Nandra et al. 1997a, b; 
Reynolds 1997; Turner et al. 1998; Iwasawa. et al. 1999; Yaqoob et al. 2002; Reynolds & Nowak 
2003; Fabian & Miniutti 2005). The bulk of the Fe K line core is likely to originate in matter 
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located much further from the central black hole than the outer disk (e.g. the outer parts of 
the optical broad line region, the putative parsec-scale torus structure, or the optical narrow 
line region). In addition to the two principal components, there are an increasing number of 
reports of cases of narrow lines at energies higher than 6.4 keV (from ionized and/or outflowing 
Fe - see Bianchi et al. 2005, and references therein) and lower than 6.4 keV (from apparently 
redshifted localized emission - Turner et al. 2002; Iwasawa, Miniutti, & Fabian 2004; Turner 
et al. 2006, and references .therein). 

Deconvolving the disk and distant-matter components of the Fe K line is extremely chal- 
lenging because the disk can also contribute to the line core emission (from radii greater than 
~ 100r 5 = 100 GM/c 2 or so). In other words, line emission from disks observed at low inclination 
angles (relative to the disk normal) is partially degenerate with line emission from distant mat- 
ter. Line emission in the wings of the profile from regions close to the black hole is smeared out 
over energy and therefore difficult to detect against the continuum. In order to derive the most 
reliable parameters and constraints for the disk and distant-matter Fe K line components, one 
must simultaneously model both, otherwise the model parameters will be biased and misleading 
(e.g. see discussions in Weaver & Reynolds 1998, Yaqoob & Padmanabhan 2004). So far only 
some of the brightest AGN have sufficient signal-to-noise ratio to produce useful constraints 
from dual-line models (e.g. see Yaqoob et al. 2002, Yaqoob & Padmanabhan 2004, and refer- 
ences therein). For the brightest sources the Fe K line core is usually unresolved even with the 
Chandra gratings because high spectral resolution and high throughput are required. Therefore, 
the location of the “distant” Fe K line emitter remains unknown. Nevertheless, decoupling the 
intensities of the broad and narrow lines (in the sense of showing that the data require both 
to be non-zero at some confidence level) would be an important step in constraining models of 
the central engine in AGN. Such decoupling has only been possible for a few of the brightest 
AGN (e.g. see Fabian & Miniutti 2005). We also note that detection of the broad Fe K line in 
type 1.5-2 AGN is even more difficult than in type 1 objects because the X-ray spectra of the 
former class are more complex and are usually affected by significant intrinsic X-ray absorption 
(e.g. see Turner et al. 1998; Weaver & Reynolds 1998; Matt et al. 2000; Guainazzi et al. 2001; 
Risaliti 2002). 

Almost all of our knowledge of the Fe K line emission in AGN comes from the Fe Ko; 
line. Although Fe K/3 line emission has been reported in a few type 1.5-2 AGN (Sambruna 
et al. 2001; Molendi, Bianchi, & Matt 2003; Matt et al. 2004), constraints on the line energy 
and equivalent width (EW) have been loose because the branching ratio, Fe K/3/Fe KM, is less 
than 0.145 (Palmeri et al. 2003, and references therein). Yet the Fe Kj3 line carries additional 
information on the ionization state of line-emitting Fe that would supplement the information 
carried by the Fe Ka line. 

In this paper we present the results of detailed X-ray CCD spectroscopy of the nearby 
(z = 0.00771, Keel 1996) Seyfert 1.9 galaxy NGC 2992 using data from Suzaku observations. 
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We show that the intensities of the broad and narrow Fe Ka lines are decoupled for the first 
time in this source. We also detected the Fe K/3 line corresponding to the Fe Ka line core and 
present measurements of its centroid energy, intensity, and EW. We describe a new method 
that utilizes the Fe K/3 line to deduce tight limits on the ionization state of Fe in the matter 
responsible for emitting the line core. For the Fe Ka; and Fe K/3 line analysis we modeled the 
0.5-10 keV spectrum self-consistently so we also report the broadband results. 

NGC 2992 has been observed by every X-ray astronomy mission since the time it was 
discovered by HEAO-1 to be one of the brightest hard X-ray AGN in the sky, with a 2-10 keV 
flux of ~ 7.2 — 8.6 x 10~ u erg cm -2 s -1 (Piccinotti et al. 1982). In more than a quarter 
of a century of X-ray observations, the hard X-ray flux has varied by over a factor of 20, 
corresponding to a range in the intrinsic 2-10 keV luminosity of 0.55 - 11.8 x 10 42 erg s" 1 
(assuming H 0 = 70 km s _1 Mpc -1 , A = 0.73, 12 = 1). During an ASCA observation in 1994, 
NGC 2992 was in its lowest continuum flux state thus far observed, with a 2-10 keV flux of 
~4 x 10 -12 erg cm -2 s~ 4 (Weaver et al. 1996). During that observation the Fe K line equivalent 
width (EW) was very high (~ 500 — 700 eV), indicating that the line intensity had not responded 
to the declining continuum. In 1997 a BeppoSAX observation found NGC 2992 to be still in 
a low continuum state, and the Fe K line EW was even higher, but then a second BeppoSAX 
observation in 1998 revealed that the source continuum had fully recovered to its bright HEAO-1 
state (see Gilli et al. 2000). In this state, the EW of the Fe K line was very low (< 100 eV) since 
the line intensity again did not respond to the change in continuum level. These measurements 
were mostly sensitive to the Fe K line core, which from the lack of variability, likely arises in a 
distant, parsec-scale reprocessor, as suggested by Weaver et al. (1996). 

The results from an XMM-Newton observation of NGC 2992 in May 2003 are unpub- 
lished. We examined those data and found that they suffered from heavy pile-up. A detailed 
analysis of the XMM-Newton data is beyond the scope of this work and will be reported else- 
where. 

NGC 2992 was the subject of an RXTE monitoring campaign consisting of 24 observa- 
tions in the period 2005 March to 2006 January. The Suzaku observations that are the subject 
of the present paper were made during the RXTE campaign and were quasi-simultaneous with 
some of the RXTE observations. The scope of the present paper is restricted to the Suzaku data 
only, but a detailed analysis of the data from the RXTE campaign will be presented elsewhere. 

The present paper is organized as follows. In section 2 we describe the observations 
and details of the data reduction. In section 3 we describe the spectral fitting analysis and 
in section 4 we discuss the detailed results of that analysis, including brief comparisons with 
historical data and results. In section 5 we describe a new method to tightly constrain the 
ionization state of Fe using the narrow Fe Ka and Fe K/3 emissions lines. In section 6 we 
summarize our results and give our conclusions. 
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2. Observations and Data Reduction 


The joint Japan/US X-ray astronomy satellite, Suzaku (Mitsuda. et al. 2006), was 
launched on 10 July, 2005. NGC 2992 was observed on three occasions in 2005, November 
and December (the observation log is given in table 1). Suzaku carries four X-ray Imaging 
Spectrometers (XIS - Koyama et al. 2006) and a collimated Hard X-ray Detector (HXD - 
Takahashi et al. 2006). Each XIS consists of four CCD detectors at the focal plane of its 
own thin-foil X-ray telescope (XRT - Serlemitsos et al. 2006), and has a field-of-view (FOV) 
of 17.8' x 17.8'. One of the XIS detectors (XIS1) is back-side illuminated (BI) and the other 
three (XIS0, XIS2, and XIS3) are front-side illuminated (FI). The bandpass of the FI detectors 
is ~ 0.4 — 12 keV and ~ 0.2 — 12 keV for the BI detector. The useful bandpass depends on 
the signal-to-noise ratio of the source since the effective area is significantly diminished at the 
extreme ends of the operational bandpasses. Although the BI CCD has higher effective area 
at low energies, the background level across the entire bandpass is higher compared to the 
FI CCDs. The spectral resolution of the XIS detectors has been continuously degrading with 
time but the current instrument response matrices employ the resolution before the onset of 
degradation, as described in Koyama et al. (2006). Although we used these standard response 
matrices for modeling the XIS data, we measured the widths of the Mn Ka lines from the 
on-board 55 Fe calibration sources (two per XIS) using the actual NGC 2992 observations in 
order to correctly interpret spectral fitting results. Details are given in section 2.3. 

2.1. HXD Data 

The HXD consists of two non-imaging instruments (the PIN and GSO - see Takahashi 
et al. 2006) with a combined bandpass of ~ 10 — 600 keV. Both of the HXD instruments are 
background-limited. We found that for the NGC 2992 observations the average on-source count 
rate in the full HXD/PIN band (~ 10 — 120 keV) was ~ 0.75 ct/s, whilst the corresponding 
estimated background count rate was 0.65 ct/s so that NGC 2992 constitutes only ~ 13% 
of the total signal. The situation is even worse for the GSO which has a smaller effective 
area than the PIN. In order to obtain reliable background-subtracted spectra, the background 
spectrum must be modeled as a function of energy and time. The background model is still 
under development and cannot yet reproduce the background spectral and temporal behavior 
with sufficient accuracy to give reliable HXD spectra for NGC 2992. Therefore, in this paper 
we will restrict our analysis to the XIS only and defer analysis of the simultaneous HXD data to 
future work when the background model systematics are lower. The observations of NGC 2992 
were in fact optimized for the HXD in terms of positioning the source at the aimpoint for the 
HXD (the so-called “HXD-nominal pointing” ) which gives a somewhat lower count-rate in the 
XIS than the “XlS-nominal” pointing but gives ~ 10% higher HXD effective area. 
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2.2. XIS Data Filtering, Cleaning, and Selection 

The XIS observations of NGC 2992 were made in so-called 3x3 and 5x5 edit modes 
which correspond to 9 and 25 pixel pulse-height arrays respectively that are telemetered to the 
ground (see Koyama et al. 2006 for details). The array is centered on the pixel containing the 
largest pulse height of an event. We combined the data from both modes as there did not appear 
to be any noticeable differences in the spectra obtained from the different modes. Data from 
version 0.7 of the pipeline reprocessing (Fujimoto et al. 2006; Mitsuda et al. 2006) were used 
for the analysis in the present paper. As is standard practice, only ASCA grades 0, 2, 3, 4, and 
6 were selected. Event files were cleaned to remove hot and flickering pixels with the FTOOL 
cleans is as updated for the XIS. Data were excluded during satellite passages through the 
South Atlantic Anomaly (SAA), and for time intervals less than 256 s after passages through 
the SAA, using the T_SAA_HXD house-keeping parameter. Data were also rejected for Earth 
elevation angles (ELV) less than 5°, Earth day-time elevation angles (DYEJSLV) less than 
20°, and values of the magnetic cut-off rigidity (COR) less than 6 GeV/c 2 . Version 0.7 of the 
reprocessing pipeline also includes CTI corrections based on available in-flight calibration and 
other refinements to the XIS energy scale. Residual uncertainties in the energy scale are on the 
order of ~ 10 eV (Koyama et al. 2006) and we confirmed this from an analysis of the onboard 
calibration line data taken during the NGC 2992 observations (see section 2.3 for details). 

The cleaning and data selection resulted in net exposure times in the range 31.6 ks 
to 41.7 ks for the three observations (see table 1). We examined XIS images made from the 
cleaned events files and found that the galaxy NGC 2993 (~ 3' from the core of NGC 2992) 
was clearly detected but its total count rate was less than 1% of the NGC 2992 count rate. 
Even though the half-power diameter of the XRT+XIS is ~ 2', NGC 2992 and NGC 2993 are 
resolved. NGC 2993 was detected by Chandra with much better spatial resolution (see Colbert 
et al. 2005). We extracted XIS lightcurves and spectra using two sets of circular regions with 
radii of 2.6’ and 4.35’ centered on NGC 2992. The smaller region excludes NGC 2993, resulting 
in negligible contamination of the NGC 2992 spectra, and the larger region corresponds to the 
size assumed by some of the standard XRT effective area files that have been released. The 
spectra from the smaller regions will be used for the principal analysis and those from the larger 
regions will be used to better estimate absolute fluxes (see below). We extracted background 
spectra using circular regions with radii of ~ 4' from source-free regions of the detectors. 

2.3. Instrument Response Files and the Mean XIS Spectrum 

For spectral analysis we used the response matrices (or “RMF” files) 
ae_xi0_20060213 . rmf , ae_xil_20060213c .rmf , ae_xi2_20060213 . rmf , and 

ae_xi.3_20060213.rmf for XIS0, XIS1, XIS2 and XIS3 respectively (binned to the 
same number of channels as the spectra). The telescope effective area files (or 
“ARF” files) used were ae_xi0_hxdnom4_20060415.arf, ae_xil_hxdnom4_20060415 . arf , 
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ae_xi2_hxdnom4_20060415 . arf and ae_xi3_hxdnom4_20060415 , arf for XISO, XISl, XIS2, 
and XIS3 respectively. The ARF files are standard files released for the HXD-nominal 
pointing position and assume a circular extraction with a radius of 2.9’. Since our 

source extraction regions had a radius of 2.6’ we derived a mean correction factor by 
comparing broadband fluxes measured from the spectra extracted using 4.35’ radii re- 
gions with those from our 2.6’ radius region. The ARF files used for the spectra from 
the larger regions were ae_xi0_hxdnom6_20060415.arf, ae_xil_hxdnora6_20060415 . arf , 
ae.jxi2_hxdnom6_20060415.arf and ae_xi3_hxdnom6_20060415 . arf for XISO, XISl, XIS2, 
and XIS3 respectively. Although the larger regions include NGC 2993, the contamination is 
small and the procedure is adequate for the present purpose. 

Some preliminary examination and spectral fitting revealed no significant evidence for 
variability in the total count rate or spectral shape within each of the three observations. 
Although there was up to ~30% variability in net count rate between observations, we found no 
statistically significant evidence for spectral variability between the observations (see table 1). 
We found that the four XIS sensors were consistent with each other to ~ 10% or better (within 
statistical errors) over the energy band 0.5 — 10 keV. However, XISO showed a systematic dip 
in the ^5 — 7 keV band (i.e. in the region that Fe K a line emission is expected) which was not 
apparent in the other three XIS sensors. Also, XISl showed a broad dip in the ~ 2 — 5 keV band 
which was not apparent in the other three XIS sensors. Although both of these effects are at 
the level of ~ 10% or less, for the sake of maximum conservatism, in the remainder of this paper 
all of the spectral analysis will refer to a single XIS spectrum combined from observations 1, 2 
and 3, using only XIS2 and XIS3 data. For analysis of the mean spectrum from XIS2+XIS3 
(integrated over the three observations), we combined the RMF and ARF files (four in total) 
into one response file with appropriate weightings, according to the exposure times of the 
data segments. The mean 0.5-10 keV background-subtracted count rate for this spectrum was 
0.4442T0.0016 ct/s per XIS, and the net exposure time was 1.08 x 10 5 s per XIS. We restricted 
the analysis to the 0.5-10 keV band in order to avoid background-subtraction systematics (in 
the 0.5-10 keV band the background constitutes ~ 3.4% of the total count rate). We obtained 
a flux correction factor for the size of the extraction region of 2.1% (as described above) for 
the mean spectrum. All fluxes and absolute line intensities quoted in this paper have been 
corrected (increased) by this factor. 

We measured the centroid energies of the Mn Ka lines from the calibration sources, 
using a background-subtracted spectrum combined from XIS2 and XIS3 data, summed over 
observations 1-3 (i.e. the same combination used for the analysis of NGC 2992), and summed 
over both calibration sources on each XIS. The expected energies of the Mn Ka \ and Mn ifa 2 
lines are 5.89875 keV and 5.88765 keV respectively (Bearden 1967). Since the Ka,\ : Ka 2 
branching ratio is 2:1, the expected centroid energy is then 5.89505 keV. Using a single Gaussian 
we measured a centroid energy of 5.906 ±0.003 keV (one-parameter, 90% confidence error). 
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Therefore, the difference between the measured and expected centroid energy is 11 ±3 eV, 
consistent with Koyama et al. (2006). Both the residual CTI correction uncertainty and the 
spectral resolution are better at the centers of the detectors (i.e. in the region of the aimpoints 
for source observations) than at the detector edges where the calibration sources are located so 
this is an over-estimate. We also measured the Mn Ka line width and obtained a FWHM of 
161.5 eV (at an energy of ~ 5.9 keV). At 6.4 keV this corresponds to a FWHM of 168.4 eV (see 
Koyama et al, 2006), or ~ 7890 km s -1 FWHM. Since the FWHM resolution at 6.4 keV assumed 
in the response matrices is 133.8 eV (from direct inspection), there is a residual FWHM width 
of 34.6 eV (or 1620 km s -1 at 6.4 keV) that is not currently accounted for in the instrument 
response function. Again, the spectral resolution is better at the center of the detectors (by 
~ 5 — 10 eV) than at the edges. 

2-4- XIS Contaminant Correction 

The XIS detectors have suffered a loss in low-energy efficiency due to absorption by a 
contaminant on the optical blocking filter. The degradation can be modeled using “effective” 
column densities of Carbon and Oxygen that can be calculated as a function of time since launch 
and position on each XIS, based on empirical fits to in-flight data obtained from non-varying 
astrophysical sources (see Koyama et al. 2006). We calculated the Carbon column density for 
each XIS and each observation of NGC 2992. For analysis of the mean spectrum described above 
we then calculated the average of the six column densities for XIS2 and XIS3, observations 1 
to 3, and obtained 3.4 x 10 18 cm -2 . Based on the current assumption that the contaminant 
is probably C 2 4 H 38 04 , the O column density was fixed at 1/6 of the value of the C column 
density. The C and O column densities were then combined with photoelectric absorption 
cross-sections from Balucinska-Church & McCammon (1992) since these were used to derive 
the empirical degradation formulae. The resulting contamination model, implemented using 
the varabs model in XSPEC, was included along with all astrophysical models of NGC 2992 
(see section 3). 

3. Spectral Fitting 

In the remainder of this paper we shall refer to analysis of the combined spectrum 
from XIS2 and XIS3 averaged over all three observations (table 1), as described in section 2.3. 
We used XSPEC vll.3.2 for spectral fitting (Arnaud 1996), utilizing data in the 0.5-10 keV 
band. The instrument response files used were described in section 2.3. The spectrum was 
binned to ~ 29 eV per bin and this resulted in > 25 ct/bin for the net spectrum over the 
0.5-9. 7 keV range so we were able to use the x 2 statistic for finding the best-fitting model and 
for the statistical error analysis. Unless otherwise stated, all statistical errors given correspond 
to 90% confidence for one interesting parameter (A^ 2 = 2.706). Galactic absorption with a 
column density of 5.26 x 10 20 cm -2 (Dickey & Lockman 1990) was included in all of the models 
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described hereafter and its inclusion will be implicitly assumed. For the Galactic absorption, 
as well as any intrinsic absorption, we used photoelectric cross-sections given by Morrison & 
McCammon (1983). The correction for the XIS contamination degradation (as described in 
section 2.4) was also included in all of the spectral fitting described below. All astrophysical 
model parameter values will be given in the rest frame of NGC 2992. 

3.1. Preliminary Spectral Fitting and the Baseline Model 

Preliminary spectral fitting confirmed the complexity in the X-ray spectrum of NGC 2992 
that has been established from observations with previous X-ray astronomy missions (see sec- 
tion 1). This complexity is illustrated in figure 1(a), which shows the ratio of the XIS data 
to a simple model consisting of a power-law intrinsic continuum and uniform intrinsic absorp- 
tion fully covering the source. There are only three free parameters in this model, namely the 
normalization of the power-law continuum, its photon index, T, and the column density, Aty. 
For this model x 2 — 1158. 2 for 321 degrees of freedom. Relative to this simple model there is 
a large excess at soft X-ray energies below ~ 2 keV and a very prominent Fe Ka emission line. 
In addition, a weaker emission line is apparent at ~ 7 keV, which could be due to Fe xxvi Ly a, 
Fe K /?, or a blend of both. From a detailed analysis described in section 5 we identify the 
weaker line with Fe K/5 and show quantitatively that contamination by any Fe xxvi Lya line 
emission is negligible. 

At least part of the soft excess in figure 1(a) is due to scattering of the intrinsic X-ray 
continuum into the line-of-sight by an extended optically-thin zone in the active nucleus (e.g. 
Turner et al. 1998). This scattered component has been imaged in NGC 2992 using ROSAT 
and Chandra (see Colbert et al. 2005, and references therein) and in some other Seyfert 1.5-2 
galaxies (e.g. Morse et al. 1995; Ogle et al. 2000; Young, Wilson, & Shopbell 2001; George 
et al. 2002; Iwasawa et al. 2003). The direct continuum along the line-of-sight is obscured so 
even if only a small fraction of the intrinsic continuum is scattered, the soft X-ray spectrum can 
easily be dominated by the scattered component. Therefore, in the next step we added another 
power-law continuum component to the model, with the photon index tied to that of the direct 
continuum and a normalization that was allowed to float. We parameterized this normalization 
in terms of the fraction, f s , of the directly observed power-law continuum normalization. In 
the optically-thin limit the Thomson depth, r es , of the scattering region (assuming spherical 
symmetry) is simply equal to f s . We also added two unresolved Gaussian components to model 
the emission lines in the Fe K band. The intrinsic widths (a^) of the two lines were tied together 
(in this case the two emission lines originate from the same atoms/ions). In the preliminary 
spectral fitting o N was fixed at a value (0.005 keV) much less than the XIS spectral resolution. 
The other parameters were the centroid energies of the two emission lines (E Fe Ka and E Fe K p), 
and their intensities (/ Fe Ka and I Fe k /3 respectively). Thus, in this more complex model there 
were a total of eight free parameters (T, Aty, f s , E Fe Ka , E Fe K/3 , I Fe Ka , / Fe K p, and a continuum 
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normalization). The best-fitting model gave x 2 = 400.6 for 316 degrees of freedom. The ratio 
of the XIS data to this best-fitting model is shown in figure 1(b). The fit is still poor and 
statistically significant residuals are still evident in the soft X-ray band below ~ 2 keV and in 
the Fe K band. Note that if we included only one Gaussian component (to model the strongest 
emission line) we obtained x 2 = 418.8 for 318 degrees of freedom, demonstrating that the weaker 
emission line is statistically significant (see section 5 for a more quantitative analysis) . We then 
allowed the intrinsic width (<t/ y) of the Gaussian emission-line components to be free and re- 
fitted the model, obtaining another significant reduction in x 2 , which dropped to 375.4 for 315 
degrees of freedom. 

Motivated by the fact that an extended soft thermal emission component has been 
detected in NGC 2992 (see Colbert et al. 2005) and in other type 1.5-2 Seyfert galaxies (e.g. 
Ptak et al. 1997; Bianchi, Guainazzi, & Chiaberge 2006, and references therein), we then added 
a component to model optically-thin thermal emission. We used the APEC model in XSPEC with 
the element abundances fixed at cosmic values (but later we allowed the Fe abundance to be free 
in order to derive constraints on it). The APEC model then involves two more free parameters, 
kT, and L A PEC , the temperature and luminosity of the plasma respectively (although we used 
the normalization of the APEC model as the actual free parameter in the fits, not its luminosity). 
In all, there were now eleven free parameters and we obtained x 2 = 354.1 for 313 degrees of 
freedom. 

There still remained some residuals in the Fe K band indicative of additional, broader, 
line emission. This might be expected from an X-ray illuminated accretion disk whose inner 
radius extends to within ~ 10 or so gravitational radii of the putative central black hole. Such 
relativistically broadened Fe K line emission has been observed in other AGN (e.g. Fabian & 
Miniutti 2005, and references therein). Indeed, it is possible that the Gaussian emission-line 
components in the above preliminary model could be modeling part of the disk line emission 
as well as the narrow core emission that might originate from more distant matter. With 
CCD data the line core can be degenerate with part of the disk line component because both 
components can have significant emission near the line rest-frame energy. Even though the 
Fe K line core may be the dominant component, it is in general necessary to simultaneously 
model both the broad and narrow components otherwise the inferred model parameters for 
one component will be biased by the presence of the unmodeled component. For example, in 
the spectral fits described above with the Gaussian widths free, the best-fitting lines may be 
broader than the true width of the Fe K line core because they may be trying to model an 
underlying broad component as well. 

Therefore, we added Fe K line emission from a relativistic disk around a (Schwarzschild) 
black hole using the diskline model in XSPEC (see Fabian et al. 1989). The diskline model 
has six parameters, namely the energy of the line in the disk frame, E Q \ the inner and outer 
radii of the disk (R- m and R out respectively, in units of gravitational radii, r g = GM/c 2 ); the 
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power-law index, q, of the radial line emissivity law ( r q ); the inclination angle of the disk normal 
relative to the observer’s line-of-sight, 9 0 b S ; and the integrated intensity of the line, /disk- The 
overall Fe K line profile in the XIS data is already dominated by the line core and we found 
that the data for any residual broad-line emission could not constrain all of the disk line model 
parameters. In particular, there is considerable degeneracy between E 0 and # 0 bs- Therefore, we 
fixed E 0 at 6.40 keV (corresponding to the rest-frame energy of Fe I Ka). In addition, we fixed 
R in and R out at 6 r g (the radius of the last stable orbit around a Schwarzschild black hole) and 
1000r ff respectively. For self-consistency we also included a disk line component corresponding 
to Fe I K/3 with no additional free parameters: R in , R ou t , q, and 9 0 b s were forced to take the 
same values as the Fe I< a line, E 0 was fixed at 7.056 keV, and the intensity of the Fe K/3 line was 
forced to be 13.5% of the Fe Ka line intensity (see section 5 for a detailed discussion on Fe K/3 
line energies and the Fe K/3/Ka branching ratio). Although there may be no clear isolated 
detection of the broad Fe K/3 line, for ionization states of Fe I to Fe xvi it will be produced , and 
excluding it from the model could affect the values of the inferred parameters of the other Fe K 
line components. Fitting the XIS data with the Fe K disk line emission components included 
in the model results in a drop in x 2 of 19.7 (to 334.4 for 310 degrees of freedom) for an extra 
three free parameters (q, 9 0 b S , and /disk)) indicating that a broad-line component is required by 
the data with a confidence level of more than 3cr. 

Finally, we added a component to model a Compton-reflection continuum using the 
model href 1 in XSPEC (see Dovciak, Karas, & Yaqoob 2004). Although only the ~ 7 — 10 keV 
band XIS data are able to constrain the reflection continuum, if any of the Fe K emission 
line components originate in Compton-thick matter, a reflection continuum will be produced. 
Excluding a reflection continuum in the model could bias the best-fitting parameters of the Fe K 
emission lines. The reflection model involves only one extra free parameter, R, which is the 
normalization of the reflection continuum relative to the case of a steady-state X-ray continuum 
illuminating a slab subtending a solid angle of 2tt at the X-ray source. The inclination angle 
between the disk normal and the observer’s line-of-sight was tied to the disk inclination angle 
for the relativistic Fe K emission-line components. Note that a disk geometry is likely to 
be inappropriate for any reflection continuum associated with Compton-thick distant-matter 
components of the Fe K line emission, but the data are not sensitive to differences in the 
reflection continuum due to geometry. When we fitted the above model to the data we found 
a best-fitting value of R — 0, and consequently no improvement in the fit. Nevertheless, we 
retained the reflection continuum in all of the spectral fitting analysis (with R free), since it 
could potentially affect the statistical error analysis for the other parameters. 

Thus, our baseline model consisted of fifteen free parameters: two normalizations (one 
for the soft, optically-thin thermal continuum, another for the intrinsic power-law continuum), 
kT, T, Nh, f s , <T/v, E Fe Ka, Ae Ka, E Fe K /3, Ae K/3, 7, #obs, Ai s k, and R. This model gave an 
excellent fit to the data (x 2 = 334.4 for 309 degrees of freedom), as can be seen from figure 2 (a) 
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(which shows the best-fitting model overlaid on the XIS counts spectrum), figure 2 (b) (which 
shows the best-fitting model and some of the individual model components), and figure 2 (c) 
(which shows the ratio of the XIS data to the best-fitting model). The best-fitting parameters 
and 90%, one-parameter statistical errors are shown in table 2. We discuss these results in 
detail in section 4. 

3. 2. Deconvolution of the Broad and Narrow Fe K Emission Lines 

Figure 3(a) shows a close-up of the Fe K region of the XIS spectrum of NGC 2992 in 
the form of the ratio of the data to a model consisting of a simple power-law continuum with 
a fully-covering absorber (i.e. figure 3(a) is simply a “zoomed’ region from figure 1(a)). Note 
that the Fe K/3 line is clearly seen here. A broad base to the prominent Fe K a line core is also 
apparent. In section 3.1 we found that when we modeled the underlying broad emission with a 
disk line, the Gaussian components that modeled the line core emission became narrower than 
for the case when the broad line emission was not modeled. In fact the best-fitting baseline 
model (table 2) indicates that the line cores of Fe Ka and Fe K/? are unresolved by the XIS data. 
In order to better illustrate the Fe K line profile due to the putative relativistic disk emission, 
we removed the disk line model components from the baseline model, froze the parameters of 
the thermal emission, reflection continuum, and Gaussian line core models and re-fitted the 
data. The resulting ratio of data to model is shown in figure 3(b) which illustrates the profile 
of the relativistic disk component of the Fe K line. 

Next, starting from the baseline model again, with its fifteen free parameters (section 3.1) 
we constructed joint confidence contours of the equivalent width (EW) of the Fe Ka line 
core versus its intrinsic width (FWHM). The 68%, 90%, 99% confidence contours are shown 
in figure 3(c) which demonstrates that the 99%, two-parameter upper limit on the FWHM 
appears to be ~ 6300 km s _1 . However, as mentioned in section 2 the response matrices 
do not include the time-dependent degradation in the XIS energy resolution. In section 2.3 
we showed that analysis of the Mn Ka calibration line data taken during the NGC 2992 
observations gave a residual FWHM width of 34.6 eV (or 1620 km s _1 at 6.4 keV) that is 
not accounted for by the current response matrices. Since the observed line width from any 
observed source is a convolution of the intrinsic line width and the instrumental broadening, 
the true 99% confidence, two-parameter upper limit of the Fe Ka line core in NGC 2992 is then 
[6300 2 — 1620 2 ]? ~ 6090 km s _1 . Note that the narrow-line widths and their one-parameter 
statistical errors in table 2 have also been corrected for this effect but the contours in figure 3(c) 
are not corrected. The relativistic disk line parameters have not been corrected. The disk line 
is then observed to be a little broader than it really is so that our lower limit on the disk 
inclination angle is slightly over-estimated and the radial emissivity parameter, q, should be 
slightly smaller. 

We also constructed joint confidence contours of the intensity of the Fe Ka disk line 
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versus the Fe K a intensity of the narrow, distant-matter line. During this process, as for fig- 
ure 3(c), all of the thirteen remaining parameters in the baseline model were free. The resulting 
contours are shown in figure 3(d), for 68%, 90%, 99%, and 3cr respectively. Remarkably, it can 
be seen that at a confidence level of greater than 3a, the disk and distant-matter Fe Kai lines 
are decoupled in the sense that the intensity of either of the two line components is non-zero 
everywhere on the 3a contour. 

We investigated the possibility that part of the disk line may be an artifact due to 
the current XIS response matrices inadequately modeling the low-energy tail of the response 
function (see Koyama et al. 2006). Again we used the Mn Ka calibration line background- 
subtracted spectrum, integrated over observations 1-3, summed over XIS2 and XIS3, using 
both calibration sources on each XIS. We modeled the main peak with a disk line and Gaussian 
line, plus a simple power-law continuum. The disk-line model parameters (see section 3.1) E 0 , 
Rjn, Rout were fixed at the values used for modeling NGC 2992 and the values of 9 ohs , and q were 
fixed at the best-fitting values shown in table 2. The disk line and Gaussian line intensities, 
as well as the peak energy and width of the Gaussian component were free parameters. By 
constructing confidence contours of the disk line intensity versus the Gaussian line intensity 
we established that the low-energy excess in the actual in-flight XIS response can only account 
for a disk line intensity that is 2% of the Gaussian line intensity, with a 99% confidence two- 
parameter upper limit of ~ 6%. Yet for the actual NGC 2992 data we obtained a best-fitting 
Fe Ka disk line intensity that was 76% of the line core (see table 2). Even the 3a confidence 
contour of the disk line intensity versus the line core intensity (see figure 3(d)) only goes down 
to a value of ~ 8% for the ratio of the disk line to line core intensity. Therefore the low-energy 
excess in the tail of the XIS response does not have a significant impact on the measurements 
of the disk line in the NGC 2992 data or its decoupling from the Fe Ka line core. 

Decoupling of the Fe Ka components from the accretion disk and distant matter has not 
been possible with previous data for NGC 2992 (see section 4.4). It is even more remarkable 
that the decoupling has been possible with CCD resolution data and for a relativistic disk line 
that is not particularly impressive in terms of signal-to-noise ratio (figure 3(b)) or EW (see 
table 2). The reason can be traced to the fact that the narrow line has excellent signal-to-noise 
and a large EW. This in turn leaves less freedom for the disk line model. While it is possible that 
there may not be two distinct physical locations of Fe Ka line emission, a single line-emitting 
region would demand a peculiar line-emissivity law as a function of radius. The strength of the 
line core compared to the broad base (which we model as a disk line) would demand that the 
emissivity fall as a function of radius and then increase again. The width of the broad base of 
the Fe Ka line (see figure 3(b)) is ~ 33,000 km s -1 FWHM, more than five times the upper limit 
on the width of the unresolved line core (figure 3(c)). Yet a uniform spherical distribution of 
matter gives line profiles which have much less contrast in the line width between the center and 
wings of the profile (e.g. see Yaqoob et al. 1993). The implications of the parameters derived 
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from the dual Fe K a line modeling will be discussed further in section 4.4 and compared with 
published results on historical data. 

4. Discussion of Results 

In this section we discuss in detail the results that were obtained by fitting the baseline 
model (section 3.1) to the XIS data (see table 2). 

4-1- Soft, Optically-thin Thermal Emission 

We obtained a temperature and luminosity of the optically-thin thermal emission com- 
ponent of kT = 0.656io!o6i keV and ~ 1.2 ± 0.4 x 10 40 erg s _1 respectively. The temperature is 
typical of that found for a similar component in other type 1.5-2 AGN (e.g. see Ptak et al. 1997; 
Bianchi et al. 2006, and references therein). Note that the errors quoted here on the temper- 
ature are statistical only and do not include systematic errors, which may be larger. Chandra 
was able to image the soft X-ray emission better than Suzaku but was still not able to spatially 
separate it from other emission components due to contamination from the PSF wings caused 
by the dominant nuclear emission (Colbert et al. 2005). Compounded with this were CCD pile- 
up issues, and lower signal-to-noise ratio of the data. Therefore Chandra obtained somewhat 
worse constraints on the temperature of the soft thermal emission. Colbert et al. (2005) give 
kT = OAllo^g keV and L = 1.4 x 10 40 erg s -1 and conclude that the bulk of the soft X-ray 
emission originates from within ~ 150 — 500 pc of the central engine. The temperature and 
luminosity of the extended soft emission component obtained by Chandra are both consistent 
with our respective values measured with Suzaku. 

Our Suzaku measurement of the temperature of the soft thermal emission in NGC 2992 
is the most reliable to date. Elvis et al. (1990) reported detecting the soft X-ray emission 
component in NGC 2992 using the Einstein HRI, but it was only detected in one out of four 
quadrants, the temperature was poorly constrained, and the luminosity was more than an 
order of magnitude greater than that obtained by Suzaku and Chandra. Although the soft 
X-ray emission imaged by the AOSMTHRI was spatially consistent with Chandra (see Colbert 
et al. 2005), the POS54THRI could not measure the temperature since it was not a spectroscopic 
instrument. The ROSAT PSPC spectral resolution was too poor to measure the temperature. 
ASCA CCD spectra had better spectral resolution but the signal-to-noise ratio of the NGC 2992 
ASCA data is lower than that of the Suzaku data (see Weaver et al. 1996). BeppoSAX could 
not constrain the soft thermal component because both the MECS and LECS instruments had 
poorer spectral resolution than the Suzaku XIS and the useful MECS bandpass did not extend 
to below ~ 0.9 keV. 
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4.2. Intrinsic Continuum and Absorption 

We obtained F = 1.57lo!o3 f° r the power-law photon index of the intrinsic continuum (see 
table 2). This is somewhat less than that obtained from BeppoSAX (r ~ 1.7 - Gilli et al. 2000) 
and Chandra (F ~ 1.8 - Colbert et al. 2005) observations. ASCA could not constrain F for the 
complex models required to explain the data (Weaver et al. 1996). Given the model-dependence 
of T and residual instrument cross-calibration uncertainties of the order of Ar ~ 0.1, we do not 
interpret the apparent differences in F amongst the different observations as necessarily having 
an astrophysical origin. 

We derived a column density of 8. 01°! x 10 21 cm~ 2 for the line-of-sight X-ray obscuration 
(see table 2) which corresponds to a Thomson depth of ~ 0.006lo!oo4 (he. the absorber is 
Compton-thin, at least in the line-of-sight). This column density is consistent with values 
obtained from previous observations with BeppoSAX (Gilli et al. 2000) and ASCA (Weaver 
et al. 1996). 

We tested the data to see if any additional Fe K absorption was required by adding a 
simple absorption-edge model (zedge in XSPEC) with the threshold energy fixed at 7.11 keV 
(corresponding to neutral Fe) but allowing the optical depth at the threshold energy to be a 
free parameter. All the other 15 free parameters in the model remained free. We obtained a 
best-fitting threshold optical depth of zero and a one-parameter, 90% confidence upper limit 
of 0.06. Since the optical depth at the Fe K edge threshold energy for the absorption already 
included in the model is ~ 0.008 (for N H =8x 10 21 cm~ 2 and the Fe abundance used), we could 
only obtain a loose upper limit on the possible over-abundance of Fe of ~ 7.5. 

Table 2 gives the observed XIS fluxes (including the effect of Galactic and intrinsic 
absorption but not the CCD contaminant) in the 0.5-2 keV and 2-10 keV bands which were 
0.16 and 1.15 x 10 -u erg cm -2 s -1 respectively. However, 'note that systematic uncertainties 
in the absolute flux calibration of X-ray detectors, as indicated by flux reference values quoted 
for the Crab Nebula observed by different instruments, are generally a few percent but can 
be as high as 10-20% in some cases (e.g. see Jahoda et al. 2006, and references therein). 
Historically, the lowest and highest 2-10 keV observed fluxes for NGC 2992 have been measured 
by ASCA in 1994 (~ 0.4 x 10 -u erg cm -2 s -1 - Weaver et al. 1996), and by HEAO-1 in 1978 
(~ 8.6 x 10~ n erg cm -2 s -1 - see Piccinotti et al. 1982) respectively. This represents a factor 
~ 21.5 in dynamic range. The Suzaku flux is at the lower end of the historically observed range, 
but it is nearly a factor ~ 3 higher than the historical minimum. 

Table 2 also gives the inferred XIS luminosities after correcting for absorption. These 
were calculated by setting the Galactic and intrinsic absorbing column densities to zero in the 
best-fitting baseline model (section 3.1). We obtained 0.7 and 1.6 x 10 42 erg s -1 in the 0.5-2 keV 
and 2-10 keV (rest-frame) bands respectively (assuming H 0 = 70 km s~ x Mpc, A = 0.73, 0 = 1). 
The historical range in 2-10 keV luminosity is therefore ~ 0.55 — 11.8 x 10 42 erg s -1 . 
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4-3. Scattered Continuum 

Table 2 shows that we measured (7.2±2.1)% for the continuum scattered into the line-of- 
sight (as a percentage of the direct line-of-sight continuum before absorption). This translates 
directly into r es = 0.072 ±0.021 for the extended, warm, electron scattering zone. The scattered 
fraction and optical depth cannot be easily compared with results from previous observations. 
The ASCA data had very poor signal-to-noise and the model in Weaver et al. 1996 did not 
include extended thermal emission. The BeppoSAX data did not have coverage of the soft 
X-ray band with sufficient resolution to measure f s (see Gilli et al. 2000). 

Colbert et al. (2005) determined from Chandra ACIS data that the scattering zone 
must be situated at more than ~ 150 pc from the nucleus since the spectrum from inside this 
region could be fitted with a power-law and absorption only. However, despite Chandra having 
much better spatial resolution than Suzaku, the Chandra data could not constrain the optical 
depth of the electron-scattering zone because the scattered continuum is spatially confused with 
the soft thermal emission and even with some of the direct central emission, due to the PSF 
wings. Spectroscopically, the Chandra ACIS data are inferior to the Suzaku data because they 
suffer from severe pile-up and have a much lower signal-to-noise ratio than the Suzaku data. 
For example, the 0.3-8 keV nuclear, and “near-nuclear” Chandra spectra analyzed by Colbert 
et al. (2005) had ~ 6,744 and ~ 14,600 counts respectively, as opposed to a total of nearly ~ 10 5 
counts (0.5-10 keV) in the XIS spectrum that we have analyzed here. Note also that pile-up in 
the Chandra data prohibited a reliable measurement of the flux and luminosity of the nuclear 
source. Thus, even though the spatial resolution of Suzaku is not as good as Chandra, the high 
signal-to-noise ratio in the Suzaku spectrum and absence of pile-up problems allowed a better 
measurement of the scattered component and hence a better estimate of the Thomson depth 
of the scattering zone. 

4-4 ■ Fe K Emission Lines 

In section 3.2 and figure 3(d) we showed that the intensities of the narrow and broad 
Fe Kct emission-line components were decoupled. The best-fitting line parameters and statistical 
errors are shown in table 2. We obtained only a lower limit on the disk inclination angle, 
#obs > 31°. The actual value is somewhat dependent on the model used for the broad line (here 
it is the simple Schwarzschild black hole disk model which does not include light bending - see 
Fabian et al. 1989). However, the overall roughly symmetric shape of the line (figure 3(b)) is a 
characteristic that excludes small inclination angles, independent of the detailed model, since 
smaller angles would give rise to an increasing asymmetry on the red side due to gravitational 
redshifting. We obtained an EW of 118 eV for the broad line, which is consistent with that 
expected from a “cold” accretion disk illuminated by the observed direct X-ray continuum (e.g. 
see George & Fabian 1991). The index of the line emissivity radial power law, q = — 1.5loi, is 
flatter than the average value for Seyfert 1 galaxies (e.g. Nandra et al. 1997a) but it is sensitive 
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to the inner and outer disk radius parameters, which we held fixed for the model-fitting (see 
section 3.1). We tried fixing q at the steeper value of —3 and obtained a lower limit on # 0 b S and 
R in of 29° and 31r g respectively. 

For the centroid energy of the Fe Ken line core (which we modeled with a Gaussian 
component), we obtained 6.407± 0.007 keV (table 2). The statistical errors on this centroid 
energy are comparable to the residual systematic errors in the absolute energy scale of the XIS 
(on the order of ~ 10 eV or so - see section 2.3, Koyama et al. 2006). For the corresponding 
Fe K/3 line we obtained a centroid energy of 7.034lo;o38 keV. The 68%, 90%, and 99% confidence 
contours of the narrow Fe K/3 line energy versus the Fe K a line energy are shown in figure 4(a), 
and figure 4(b) shows corresponding contours for the line intensities of Fe Ka and Fe K/3. 
Note that both sets of contours were constructed from the baseline model (section 3.1) allowing, 
in each case, all thirteen of the remaining parameters to be free. The measured energies and 
statistical errors on the Fe Ka and Fe K/3 line energies are a function of calibration uncertainties, 
Doppler and gravitational shifts, and the ionization state of Fe. Further analysis utilizing 
information from the Fe K/3 line, including the Fe K/3/Ka ratio, is described in section 5. 

The measured intensity of the narrow (distant-matter) Fe Ka line, In = 2-5-qa x 
10~ 5 photons cm -2 s -1 (see table 2), is a factor of ~ 20 times larger than the line intensity ex- 
pected from a fully-covering Compton-thin spherical shell with the cosmic abundance of Fe used 
in this paper, the measured line-of-sight column density (8.0lo| x 10 21 cm -2 ), and the measured 
intrinsic 2-10 keV luminosity (see table 2). Using the expressions in Yaqoob et al. 2001 (see 
also Krolik & Kallman 1987), we obtained a predicted intensity of 1.1 x 10~ 6 photons cm -2 s _1 . 
Therefore, if the line emitter is Compton-thin, one or more of the following must be true: (i) 
the Fe abundance is higher than the assumed cosmic value, (ii) the column density out of the 
line-of-sight is greater than that measured along the line-of-sight, and (iii) there are time delays 
between changes in the amplitude of the X-ray continuum luminosity and the response of the 
narrow Fe Ka line so that the observed line intensity may correspond to a time-averaged con- 
tinuum illumination. On the other hand, the EW of 163126 e V measured for the narrow Fe Ka 
line is also compatible with that expected from a line observed from the X-ray illuminated inner 
surface of a Compton-thick torus-like structure (although the historically averaged illuminating 
luminosity may need to be higher than the luminosity directly observed, e.g. see Makishima 
1986; Ghisellini, Haardt, & Matt 1994). However, if the Fe Ka line core were produced in 
Compton-thick distant matter we would expect to observe the so-called “Compton shoulder” 
on the red side of the line core. The Compton shoulder would have manifested itself as an 
asymmetry in the E Fe k /3 versus E Fe k a contours shown in figure 4(a), in the sense that they 
would be more elongated towards the lower energies (e.g. see Yaqoob et al. 2005 and references 
therein). Also, upper limits on the expected Compton-reflection continuum appear to be too 
low (section 4.5). Thus, the distant-matter Fe Ka line is likely to originate in Compton-thin 
matter. 
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Table 2 and figure 3(c) show that the Fe Ka line core (narrow component) is unresolved 
by the XIS, and we obtained a 90% confidence, one-parameter upper limit of 4630 km s _1 on 
the FWHM (or < 42 eV). These upper limits on the line width have been corrected (as 
described in section 3.2) for the fact that the current instrument response matrices do not 
model the time-dependent degradation of the XIS spectral resolution (see Koyama et al. 2006). 
The measured upper limits on the line width can be compared with a theoretical estimate of 
FWHM ~ 760^(M & /r pc ) km s -1 (assuming a virial relation and an r.m.s. velocity dispersion of 
v / 31 / fwhm/2; e.g. see Netzer 1990). Here, M 8 is the central black-hole mass in units of 1O 8 M 0 , 
and r pc is the distance (in parsecs), of the line-emitting structure from the central mass. The 
mass of the central black hole in NGC 2992 has been estimated to be 5.2 x 10 7 M 0 from stellar 
velocity dispersions (Woo & Urry 2002, and references therein). Therefore, the measured upper 
limit on the FWHM corresponds to r pc > 0.034, or r > 40.2 light days. For comparison, broad 
Pa/5 line emission hats been observed in NGC 2992 with a FWHM of ~ 1200 km s _1 (Rix 
et al. 1990) indicating ~ 250 light days for the location of the putative hidden BLR (see also 
Veilleux, Goodrich, & Hill 1997). 

It is pertinent to ask whether the narrow Fe Kcr line core has varied in intensity compared 
to historical data (see Weaver et al. 1996; Gilli et al. 2000, and references therein). Given the 
large historical amplitude of variability of the continuum luminosity (greater than a factor of 
20), the response of the narrow Fe K a line to the continuum would give some clue about the 
size of the line-emitting structure. In order to answer the line variability question one needs 
to fit spectra with broad and narrow Fe Ka' line models simultaneously, but this has not been 
done in the literature. We have fitted historical ASCA and BeppoSAX data with the dual Fe K 
line model used for the Suzaku data and have found that at the 99% confidence level, the ASCA 
(1994) and BeppoSAX (1997) data are consistent with no variability of the narrow Fe Ka line, 
and the BeppoSAX (1998) data cannot rule out variability up to a factor of ~ 5. This is a 
factor of ~ 4 less than the continuum variability and supports an origin of the narrow Fe Ka 
line core in distant matter. At 99% confidence, the intensity of the broad (disk) Fe Ka line 
during the low-state ASCA and BeppoSAX (1997) observations is consistent with no variability 
compared to the Suzaku measurement, but the high state BeppoSAX (1998) data cannot rule 
out a factor of up to ~ 15 variability. Therefore the relation between the broad line intensity 
and continuum luminosity remains open to future investigation. We note that in the 1997 
BeppoSAX data (when NGC 2992 was in a near-minimum flux state), there were additional 
residuals at ~6.8 keV (rest-frame) relative to the dual-line model possibly due to emission from 
Fe XXV and/or Fe XXVI (see discussion in Gilli et al. 2000). 

4-5. Compton- Reflection Continuum 

If the broad Fe K line originates in an X-ray illuminated accretion disk then basic physics 
predicts a specific Compton-reflection continuum shape and amplitude commensurate with the 
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Fe K line parameters (e.g. George & Fabian 1991; Magdziarz & Zdziarski 1995) and the disk 
ionization state. We used a simple disk reflection model (see section 3.1) to parameterize the 
total reflection continuum from all possible locations (e.g. the accretion disk and any Compton- 
thick distant matter) because only a small portion of the XIS bandpass (~ 7 — 10 keV) can 
possibly constrain it. The best-fitting value of the so-called reflection fraction, R, was 0. R 
is the normalization of the reflection continuum relative to that expected from a steady-state 
X-ray illuminated neutral disk subtending a solid angle of 2tt at the X-ray source. Although 
the XIS data cannot constrain the reflection continuum very well, we included it in the model 
because omitting it could potentially affect the inferred values of some of the other model 
parameters, especially those of the broad disk line. We found that the upper limits on R were 
strongly dependent on the disk inclination angle, principally because of the angle-dependence 
of the relativistic line profile shape (the angular dependence of the reflection continuum shape 
in the restricted energy band is negligible). Therefore, we constructed joint 68%, 90%, and 99% 
confidence contours of R versus 9 0 b S , and these are shown in figure 5. It can be seen that the 
90% confidence, two-parameter upper limit on R varies from ~ 0.6 to 0.9 as 6 0 b S varies from its 
lower limit of ~ 30° to ~ 80° respectively (the spectral fit is not well-behaved as we approach 
~ 90°). For comparison with values in the literature, the 90%, one-parameter upper limits on R 
are 0.60 and 0.78 for fixed inclination angles of 9 obs — 43° (the best-fitting value) and 0 O bs = 60° 
respectively (see also table 2). 

We checked that the limits on the Compton-reflection continuum were not sensitive to 
background subtraction systematics by constructing a completely new background spectrum 
from rectangular regions from off-source locations in XIS2 and XIS3 and summed over the 
three observations. Note that the nominal background spectrum used thus far was constructed 
from circular regions centered at different positions and as far as possible relative to the new 
background regions (see section 2.2). With the new background spectrum in place, the baseline 
model with fifteen free parameters was re-fitted and there was no statistically distinguishable 
change in any of the parameters. The one-parameter, 90% confidence upper limit on R for 
0 ohs = 43° remained at 0.60 and the corresponding upper limit for 9 0 b s = 60° decreased from 
0.78 to 0.76. Note in figure 5 how the data strongly reject models with small inclination angles 
($ 0 bs < 20° is ruled out at > 99% confidence). This is further confirmation that the narrow and 
broad Fe K lines are decoupled because in poorer quality data a low-inclination angle disk line 
model can fit the narrow core of the Fe K line profile even if the line core has a non-disk origin. 

We also checked that the disk line intensity was not sensitive to the amplitude of the 
reflection continuum in the fit by constructing confidence contours of the disk line intensity, 
Fdisk) versus R (and reverting to the nominal background spectrum for consistency with the 
earlier analysis). We found that that the 99% confidence lower limit on 7 disk was fairly flat and 
non-zero over the range 0 < R < 3, and was never less than 0.4 x 10“ 5 photons cm -2 s -1 . 

Under the stated assumptions of our model we would expect R= 1. The one-parameter 
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upper limits on R discussed above are not necessarily inconsistent with the measured EW 
(1181H eV) of the broad Fe K line. The small bandpass leverage may bias R, and our model 
for the reflection continuum does not include relativistic smearing (which would tend to lower 
the effective value of R). The model also assumes a cold, neutral disk with a cosmic Fe 
abundance. However, our upper limits on R do imply that there is unlikely to be a significant 
Compton-reflection continuum from the structure that produces the core of the Fe K fluorescent 
line emission (e.g. the putative obscuring torus), which in turn implies that the structure is 
Compton-thin. 

Unfortunately the PDS data from both of the BeppoSAX observations, including the 
high-state observation in 1998, were too noisy to obtain useful measurements or limits on a 
Compton-reflection continuum (see Gilli et al. 2000). Although Colbert et al. (2005) claimed to 
measure a Compton-reflection continuum in the off-nuclear region of NGC 2992 using Chandra , 
the data only extended up to 8 keV, the signal-to-noise was poor, and the data suffered from 
severe pile-up. The models fitted were complex in relation to the quality of the data and the 
uniqueness of the deconvolution of the different spectral components under such circumstances 
is an issue that was not addressed. In the future we expect the background modeling and sub- 
traction S 3 ^stematics for the HXD onboard Suzaku to improve and we will revisit the challenge 
of measuring the Compton-reflection continuum in NGC 2992. Nevertheless, it is important to 
note that the detection of the relativistic disk line in the XIS data is insensitive to uncertainties 
in the reflection continuum: there was no preferred model solution in which a non-zero value of 
R forced the relativistic disk line intensity to zero. 

f.6. Upper Limits on the Ni/Fe Abundance 

We investigated Ni Ka line emission in the XIS data by adding another Gaussian line 
component to the baseline model (section 3.1). The intrinsic width of the additional Gaussian 
was tied to that of the Fe Ka and Fe Kfl lines. Since visual inspection revealed little evidence 
for Ni Ka line emission we fixed the center energy of the Gaussian at the expected (rest-frame) 
energy of 7.472 keV (Bearden 1967). Thus, only the intensity of the Ni Ka line was a free 
parameter, giving a total of sixteen free parameters in the model. The results are shown in 
table 2, and it can be seen that a statistically significant detection was not obtained. The 
decrease in the fit statistic was only A% 2 ~ 1.1. We obtained a 90% confidence, one-parameter 
upper limit on the EW of 28 eV. 

If the Fe Ka and Ni Ka lines are formed in Compton-thin matter, the ratio of their inten- 
sities depends principally on the respective ratios of the fluorescence yields, the Kfl/Ka branch- 
ing ratios, the photoelectric absorption cross-sections, and the Fe and Ni abundances. Using 
Bambynek et al. (1972) for the fluorescence yields and Kfl/Ka branching ratios, and Henke 1 
for photoelectric absorption cross-sections, we evaluated equation (3) in Yaqoob et al. (2001) for 

1 www-cxro.lbl.gov 
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Fe and Ni and obtained the theoretical ratio of line intensities: I N ; Ko ,/I Fe k q = 0.96 (A Ni /A Fe ). 
From the XIS data we constructed a joint 90% confidence contour of In; kq versus I Fe Ka and 
obtained an upper limit on I N ; Ka/Ipe k q of 0.176. Thus we obtained a 90% confidence upper 
limit on the Ni/Fe abundance ratio of 0.183. For comparison, the Ni/Fe abundance ratio from 
the Anders & Grevesse (1989) solar values is 0.038 (i.e. a factor 4.8 smaller than our upper 
limit). If the Ka lines were observed in reflection from Compton-thick matter then the expected 
I]Ni Ka / Ipe Ka ratio would be less than the Compton-thin case (due to absorption of the Ni Ka 
by Fe) and our upper limit on the Ni/Fe abundance would be larger. 

5. Precision X-ray Spectroscopy with the Fe Ka and Fe K/3 Lines 

In this section we will constrain the ionization state of Fe responsible for the narrow Fe K 
line emission in a robust manner by using as much information as possible from the spectrum 
and by utilizing the latest atomic physics data and associated uncertainties. Firstly, we can 
make an immediate deduction from the very fact that the Fe K/3 fluorescence line is detected 
at all. That is, Fe must be less ionized than Fe xvn in the line-emitting matter responsible for 
the narrow Fe K line emission because Fe xvil and ions with higher ionization states have no 
M shell electrons. We can constrain the ionization state further by taking advantage of the fact 
that we measured the Fe Ka and Fe K/3 line energies with sufficiently small statistical errors. 
We cannot do this directly from the measured line energies (e.g. using figure 4(a)) because they 
could be affected by residual XIS energy scale uncertainties, possible gravitational redshifts, 
and possible Doppler shifts. We will describe a method which enables us to place an upper 
limit on the dominant ionization state of Fe that is independent of the residual uncertainty in 
XIS gain , and does not require any knowledge of the possible gravitational and Doppler shifts. 
Our result will depend only on the uncertainty in the absolute offset in the XIS energy scale, on 
which we can confidently place robust limits in general (Koyama et al. 2006), and for for every 
observation using the onboard calibration sources. However, we first justify our confidence that 
the Fe K/3 line suffers from negligible contamination from possible Fe I<a line emission. 

In figure 4(a) we showed the 68%, 90%, and 99% confidence contours of E Fe k /3 versus 
E Fe Ka measured from the XIS data. The fact that the contours are fairly symmetric about the 
best-fit and are not elongated towards lower Fe K/3 energies implies that there is little or no 
blending from Fe xxvi Lya, whose rest-frame centroid energy would be at 6.966 keV (e.g. Pike 
et al. 1996). We added another Gaussian line component to the baseline model (section 3.1) 
at a fixed energy of 6.966 keV, with its intrinsic width tied to that of the Fe Ka and Fe K/3 
lines and obtained a best-fit Fe xxvi Lya- EW of 0 and a 90%, one-parameter upper limit of 
18.5 eV. 

The 68%, 90%, and 99% confidence contours of I Fe K /? versus I Fe k q measured from the 
XIS data are shown in figure 4(b). The lower and upper dashed lines correspond to 7p e Kp/d Fe Ka 
ratios of 0.12 and 0.135 respectively. Palmeri et al. (2003) give a detailed discussion of the theo- 
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retical and experimental values of the Fe K/3/Fe Ka branching ratio. They compare theoretical 
values of the ratio in the literature for Fe I to Fe IX calculated using different codes and as- 
sumptions and they compare these theoretical values with several experimental measurements 
for Fe I. Palmeri et al. (2003) show that for Fe I, except for one experimental value, the re- 
mainder of the theoretical and experimental values of the Fe K/3/Fe Ka branching ratio lie in 
the range ~ 0.12 — 0.135. Including the highest experimental value, the range is ~ 0.12 — 0.145 
corresponding to a dispersion of ~ 20%. For higher ionization states of Fe the branching ratio 
increases, and Palmeri et al. (2003) show that the theoretical value for Fe IX can be as high as 
~ 0.17. The actual value of the observed ratio of the Fe K/3 to Fe Ka line depends a little on 
the optical depth of the medium in which the lines originate since the two lines have a slightly 
different albedo. 

New calculations for the energies of fluorescent Fe Ka and Fe K/3 transitions, fluores- 
cent yields, and rates, for all charge states of Fe have been presented by Palmeri et al. (2003), 
Mendoza et al. (2004), Bautista et al. (2003), and Kallman et al. (2004). These works also 
include a comparison with previous calculations (and measurements, where appropriate) in the 
literature. In figure 4(a), for Fe i-Fe ix, we have plotted the expected Fe K/3 line energies 
against the expected Fe Ka line energies (solid, filled circles), overlaid on the measured Fe K/3 
line energy versus Fe Ka line energy confidence contours. The energies were calculated from 
Palmeri et al. (2003), weighted according to the appropriate branching ratios for Ka i and Ka 2 , 
and for Kfii and Kf3 2 where necessary. As discussed in detail by Kallman et al. (2004), and 
references therein, for ionization states higher than Fe ix the number of transitions increases, 
and a single line centroid energy for comparison with observations becomes less and less mean- 
ingful since the fluorescent line profile shape then depends on the details of the excitation 
mechanisms in play. One must then compute the detailed photoionization equilibrium in order 
to derive theoretical line profiles to compare with observations. From a comparison of results 
of theoretical calculations with laboratory measurements (solid state for Fe I, EBIT results for 
Fe ix), Kallman et al. (2004) assess the systematic error in the Fe K fluorescent line energies 
to be ~ 0.5 mA for Fe I and ~ 2 mA for Fe ix. We note that in eV, a systematic error of AA 
in mA, is AE(eV) ~ 3.3AA and AE'(eV) ~ 4.0AA at 6.4 keV and 7.0 keV respectively. 

It is important to note from the filled circles in figure 4(a) that although the theoretical 
Fe K/3 line energy increases along the sequence Fe I to Fe ix, the Fe Ka line energy does not 
increase for these ionization states. This was pointed out by Palmeri et al. (2003) and Kallman 
et al. (2004) and the results are essentially a refinement of the results of Makishima (1986) 
which are universally quoted to describe a monotonic increase in the Fe Ka line energy with 
increasing charge state of Fe. However, Makishima (1986) never gave numerical values for the 
Fe Ka line energies, and for the ionization states Fe I to Fe IX it is not possible to extract 
the line energy values from their plot with sufficient precision to compare with the values from 
Palmeri et al. (2003). Although the decrease in the Fe Ka line energy in going from Fe I to Fe ix 
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is less than 10 eV, the XIS data for NGC 2992 are sensitive to these small energy differences. 
The statistical and systematic errors on the Fe K a line energy measured from the NGC 2992 
XIS data are comparable to the difference in the theoretical values of the Fe I and Fe IX line 
energies. The values of the Fe Ka and Fe K/3 line energies, the associated statistical errors, and 
confidence contours measured for NGC 2992 are not only a function of the predominant charge 
state of Fe and residual calibration uncertainties, but are also affected by possible Doppler and 
gravitational shifts. In order to make inferences on the ionization state of Fe from the XIS data 
we must consider all of these factors, including uncertainties in the theoretical line energies, so 
we cannot simply deduce the ionization states allowed by the data from figure 4(a). 

In order to constrain the dominant ionization stage of Fe further, we can recast the 
E Fe k/3 versus E Fe K a contours shown in figure 4(a). For a given ionization stage of Fe there are 
a pair of “true” (or theoretical) values of E Fe K « and E Fe k/ 3 - Suppose that the observed and 
true values of E Fe Ka and E Fe k /3 are linearly related by 

Ei- Fe Ka (observed) = A + Bx E i>Fe Ka (true) (1) 

and 

E k Fe k/3 (observed) = A + B x E i>Fe K /3 (true), (2) 

where i = 0 to 25 represents the ionization state of Fe (i.e. the number of electrons missing 
in the ion relative to the neutral atom), A represents any residual calibration uncertainties in 
the offset of the absolute energy scale and B includes contributions from residual uncertainties 
in the XIS gain as well as possible gravitational and/or Doppler shifts. Certainly, over the 
restricted energy range being considered here, any higher order terms in the gain uncertainty 
are negligible (see Koyama et al. 2006), and the gravitational and Doppler terms only make 
a linear contribution. Therefore, the linear equations above are an excellent representation of 
the relations between the true rest-franre line energies and the final line energies measured by 
the XIS. 

Since A and B are the same for the Fe Ka and Fe K/3 lines, we can solve for A and B 
for a given ionization stage of Fe for any given pair of measurements of Il/pe k q (observed) and 
E it Fe kp (observed). Thus we can recast a contour of the measured energies E Fe k/3 and E Fe Ka 
for a given confidence level into a contour of B versus A for every i (i.e. we can generate up to 
25 contours from figure 4(a)). Since we have independent information on the calibration of the 
absolute energy scale of the XIS data, we can then rule out any ionization states that yield B 
versus A contours that give solutions for A that are inconsistent with the known accuracy of 
the XIS energy scale. 

The 99% confidence contours of B versus A calculated as described above are shown in 
figure 6 for Fe I to Fe IX, using the expected Fe Ka- and Fe K/3 line energies calculated from 
Palmeri et al. (2003). These energies were also shown as filled circles in figure 4(a). The vertical 
dashed lines in figure 6 represent the uncertainty in the absolute energy scale offset of the XIS 
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(i.e. they enclose the possible range of values of A). Although we have stated that this is 
nominally ~ ±10 eV, we will base our conclusions on the very conservative assumption that \A\ 
is less than 50 eV. The vertical dashed lines are therefore at A = ±50 eV. It can be seen from 
figure 6 that the contours for Fe VII to Fe IX are ruled out regardless of the value of B. We can 
go further and take into account the uncertainties in the theoretical values of the Fe K a and 
Fe K/3 line energies. The B versus A contours depend on E Fe K p—. E Fe Ka and we considered the 
two extreme cases, corresponding to the theoretical difference in these energies hypothetically 
being in error by —2mA and ±2mA. We recalculated the B versus A confidence contours for 
each of these cases for Fe I to Fe IX. The cases for a theoretical error of —2mA and ±2mA gave 
contours which ruled out ionization states higher than Fe v and Fe vn respectively. Therefore, 
we can say that at 99% confidence the narrow core of the Fe Ka line in NGC 2992 originates 
in matter in which the predominant ionization state of Fe lies in the range Fe I— Fe VII. We 
emphasize that this conclusion requires no knowledge of the uncertainty in the XIS gain, or of 
any possible gravitational or Doppler shifts experienced by the photons leaving the source. 

6. Summary and Conclusions 

We have presented detailed X-ray spectroscopy in the 0.5-10 keV band of the Seyfert 1.9 
galaxy NGC 2992 with the Suzaku XIS. Here we summarize the results. 

1. The mean observed 2-10 keV flux of NGC 2992 during the three Suzaku observations in 
2005 November and December was ~ 1.2 x 10 -11 erg cm -2 s _1 , a factor of ~ 3 higher than 
the historic minimum flux, and a factor of ~ 6 lower than the historic maximum. The 
absorption-corrected intrinsic 2-10 keV luminosity was 1.6 x 10 42 erg s _1 . 

2. We detected the extended, optically-thin soft X-ray thermal emission found with previous 
X-ray astronomy missions, and obtained a temperature of kT = 0.656io.o6i keV. This is 
the most precise measurement of kT for NGC 2992 to date although the errors are purely 
statistical and do not include any possible systematic effects. We obtained a luminosity 
for the soft emission component of ~ 1.2 ±0.4 x 10 40 erg s _1 which is consistent with that 
deduced from Chandra data (Colbert et al. 2005). 

3. We modeled the intrinsic X-ray continuum with a power law and obtained F = 1. 571^03 
for the power-law photon index. 

4. We derived a column density of 8.0 Aoil x 10 21 cm" 2 for the line-of-sight X-ray obscuration, 
which corresponds to a Thomson depth of 0.006±8;88| (i.e. the absorber is Compton-thin, 
at least in the line-of-sight). This column density is consistent with values obtained from 
previous BeppoSAX and ASCA observations. 

5. We measured r es = 0.072 ± 0.021 for the optical depth of the extended, warm, electron 
scattering zone which is responsible for scattering (7.2 ±2.1)% of the intrinsic hard X-ray 
continuum into the line-of-sight. 
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6. We detected broad and narrow Fe K line emission. When we modeled the broad and 
narrow lines with line emission from a relativistic disk around a Schwarzschild black hole 
and a Gaussian profile respectively, we found that the intensities of the broad and narrow 
lines were decoupled at a confidence level of greater than 3cr. We obtained an EW of 
llSlgi eV, a radial line-emissivity law of r - 1 - 5 (-°- 8 .+ 5 - 1 ) ; and a lower limit of 31° on the 
disk inclination angle from the disk-line model. The actual values are model-dependent, 
however. 

7. The narrow Fe K line core EW was eV. The line core was unresolved by the 

XIS and we obtained a 90% confidence, one-parameter upper limit on the FWFIM of 
4850 km s -1 . However, as discussed in section 2.3 the current XIS response does not 
model the degradation in spectral resolution. After correction for this effect, as described 
in section 3.2, we get an upper limit of 4630 km s -1 (or a ^ < 42 eV). 

8. The absolute flux in the narrow line is a factor of more than 20 larger than that expected 
from a fully-covering spherical shell illuminated by the measured intrinsic 2-10 keV lu- 
minosity, having a cosmic abundance of Fe and a uniform column density equal to that 
measured in the line-of-sight. Therefore, roughly speaking, the product of the solid angle 
subtended by the line-emitting structure at the source, its angle-averaged column density, 
its Fe abundance, and the mean historical 2-10 keV luminosity averaged over timescales 
of the order of the structure light-crossing time, must be > 20. 

9. We detected the Fe K/3 fluorescent line corresponding to the Fe Ka line core. Taking 
into account residual uncertainties in the XIS energy scale calibration and conservative 
uncertainties on the theoretical energies of the Fe K a and Fe K/3 transitions for different 
ionic species of Fe, we were able to constrain the predominant ionization state of Fe 
responsible for producing the narrow Fe K line core to be Fe I-Fe VII. Ionization states of 
Fe VIII or higher are ruled out at 99% confidence. This is the most precise determination 
to date of the ionization state of the matter emitting the Fe K line core and was achieved 
using a method that requires no knowledge of any uncertainties in the XIS gain and any 
possible gravitational and Doppler shifts. 

10. We were able to decouple the disk and distant-matter Fe Ka emission lines because the 
continuum flux was low, giving a large EW for the line core, with good signal-to-noise 
ratio. The large core EW also enabled the Fe K/3 line to be measured with a good signal- 
to-noise ratio, giving redundant information on the width of the distant-matter lines. 
This demonstrates the importance of measuring the distant-matter Fe K line emission 
parameters in order to measure those of the relativistic Fe K line. 

11. Although the XIS data do not have the bandpass to measure the Compton-reflection 
continuum associated with the broad Fe K line, we included it in all of the spectral analysis 
with a parameterized normalization, R, relative to that expected from a steady-state 
illumination of a disk subtending a solid angle of 27r at the X-ray source. We obtained 
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upper limits on R as a function of the disk inclination angle (constrained by the relativistic 
Fe K line). We found (at 90% confidence for two parameters), R < 0.9 for an inclination 
angle of 80°, going down to R < 0.6 when the inclination angle is at its lower limit of 31°. 
The fact that these upper limits on R are less than 1 implies that the narrow Fe K line 
core is likely to originate in Compton-thin matter. 

In the future we will refine our modeling by utilizing the Suzaku HXD data when the 
background model has achieved the goal of 1% systematic uncertainties. The HXD data will 
better constrain the Compton-reflection continuum components. We will also improve model 
constraints by including data from XISO and XISl when the calibration has advanced further. 
Future observations with higher spectral resolution, such as gratings or calorimeters, will be 
able to resolve the narrow core of the Fe K a line and thereby better constrain the location of the 
line-emitter. In addition, improved throughput will enable better measurement of the profile 
of the broad Fe I< line. Further monitoring of NGC 2992 to measure variability in the Fe K 
line components and in the reflection continuum as the continuum undergoes large excursions 
in amplitude promise to reveal important clues on the accretion process and on the structure 
of the central engine beyond the accretion disk. 

We thank the Suzaku team and operations staff for their hard work and dedication 
that has made these observations of NGC 2992 possible and that has facilitated the cal- 
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Table 1. NGC 2992 Suzaku Observation Log 


Obs 

Start 0 


End“ 


Count Rate 6 

Exposure 


(UT) 


(UT) 


(ct/s/XIS) 

Time 6 (s) 

1 

6/11/2005 

14:16:51 

7/11/2005 

14:01:50 

0.3888 ± 0.0027 

34664.5 

2 

19/11/2005 

21:42:59 

20/11/2005 

23:24:45 

0.5096 ± 0.0032 

31641.5 

3 

13/12/2005 

10:15:37 

14/12/2005 

12:08:38 

0.4408 ± 0.0026 

41681.0 


° Start and end times correspond to the time tags of the first and last photons respectively in the cleaned and 
filtered events files combined from XIS2 and XIS3. 6 Mean count rates and exposure times are for XIS2 and 
XIS3 per XIS, in the 0.5-10 keV band. 
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Fig. 1. (a) Ratio of the XIS2+XIS3 spectrum (integrated over observations 1 to 3) to a model consisting 
of Galactic absorption, an intrinsic power law continuum, and a uniform intrinsic absorber, (b) As (a), 
but now the model has additional components corresponding to an unabsorbed power-law component (for 
example due to scattering), and narrow (unresolved) Gaussian Fe I(a and Fe K j3 emission-line components. 
See section 3.1 for details. 


31 




0.5 1 2 5 10 

Observed Energy (keV) 


Fig. 2. (a) The XIS2+XIS3 counts spectrum (integrated over observations 1 to 3) overlaid with the 

best-fitting baseline model (solid line). The inset shows a close-up of the region containing the Fe Kq and 
Fe K/3 emission lines. The C/O XIS contamination is included in the model, (b) The best-fitting baseline 
model (see section 3.1 and table 2). The soft X-ray optically-thin thermal continuum and the broad Fe Ka 
and Fe K/3 disk-line components are also shown separately. The C/O XIS contamination model has been 
removed in the plot but Galactic absorption is included, (c) The ratio of the data in (a) to the best-fitting 
baseline model. 
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Fig. 3. (a) A close-up of the XIS spectral region in figure 1(a) containing the Fe Ka and Fe K/3 lines. 
The vertical dotted lines correspond to the observed-frame energies of Fe I Ka and K/3 at 6.400 keV 
and 7.056 keV in the rest-frame respectively. All of the spectral analysis in the present paper refers to a 
mean spectrum constructed from XIS2+XIS3, integrated over observations 1 to 3. (b) The ratio of the XIS 
spectral data to the baseline model re-fitted after removing the relativistic Fe Ka and Fe K(3 disk-line com- 
ponents, whilst freezing the optically-thin thermal continuum, narrow Fe I< line, and Compton-reflection 
continuum parameters at their best-fitting values (see section 3.2 for details), (c) The 68%, 90%, 99% 
confidence contours of the equivalent width (EW) versus the intrinsic width (FWHM) of the narrow, dis- 
tant-matter Fe Ka line based on the best-fitting baseline model (section 3.1). The contours have not been 
corrected for the degradation in XIS spectral resolution (see section 3.2). (d) The 68%, 90%, 99%, and 3a 
confidence contours of the Fe Ka relativistic disk-line intensity (/disk) versus the intensity of the narrow, 
distant-matter Fe Ka line (/narrow = /tv) based on the best-fitting baseline model (section 3.1). It can be 
seen that even at a confidence level of 3a, the disk and distant-matter Fe Ka lines are decoupled in the 
sense that the intensity of either of the two line components is non-zero everywhere on the 3a contour. 
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Fig. 4. (a) The 68%, 90%, and 99% confidence contours of Ep e k /3 versus Ep e K«(in the NGC 2992 

rest-frame), measured from the mean XIS spectrum (see section 4.4). The filled circles are theoretical 
rest-frame energies of the Fe K/3 and Fe Ka emission lines from Palmeri et al. (2003). (b) The 68%, 90%, 
and 99% confidence contours of Jp e k/3 versus Ip e Ka(in the NGC 2992 rest-frame), measured from the 
mean XIS spectrum (see section 5). The lower and upper solid lines correspond to Fe K/3/Fe Ka branching 
ratios of 0.12 and 0.135 respectively (see section 5). 
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Fig. 5. The 68%, 90%, 99% confidence contours of the relative Compton-reflection fraction (R) versus the 
inclination angle (0 o bs) of the disk normal to the observer’s line-of-sight, based on the best-fitting baseline 
model (section 3.1). A value of R = 1 corresponds to a, Compton-reflection continuum normalization 
expected from a disk subtending a solid angle of 2ir at a central non- varying X-ray source. The inclination 
angle parameters of the disk for the Compton-reflection continuum and relativistic disk Fe Ka line model 
components were tied together. See section 4.5 for details. 
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Fig. 6. The 99% confidence contours of B versus A calculated from the 99% contour of Ep e versus 
Fpe Ka in figure 4(a) (which was derived directly from the XIS spectrum of NGC 2992), for Fe I to Fe IX. 
The method of transformation of the E-p e K 0 versus Ka contour is described in section 5. The Fe Ka 
and Fe Kj3 fluorescent line energies were calculated using atomic data in Palmeri et al. (2003) and are 
also shown as filled circles in figure 4(a). A represents residual uncertainties in the offset of the absolute 
energy scale which is believed to be on the order of ~ 10 eV, but the dashed vertical lines are placed very 
conservatively at ±50 eV. B represents the product of any gravitational shift, Doppler shift, or residual 
uncorrected gain factors. Contours of B versus A that do not have any portion lying inside A = ±50 eV 
are ruled out at 99% confidence. 
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Table 2. XIS Spectral Fitting Results for NGC 2992 


X 2 / degrees of freedom 

334.4/309 

kT (keV) 

0.658±g;®|S 

T'Apec (10 40 erg s -1 ) 

i.i8 ±8:11 

r 

l.ouy_ 0 027 

N h (10 21 cm- 2 ) 

7.99i°:f 5 

f s (scattered fraction) 

0.073±g;^i 

9 0bs (degrees) 

>31 

q 

-1 >$+ 5 - 1 

/disk [Fe Ka] (10~ 5 photons cm -2 s _1 ) 

i-9±?;S 

EWdisk [Fe Ka] (eV) 

llSlei 

En [Fe Ka] (keV) 

6 407 + ° 007 
D '^ u ‘ -0.007 

o N (keV) a 

< 0.042 

FWHM (km s- 1 ) a 

< 4630 

/n [Fe Ka] (10~° photons cm -2 s _1 ) 

2-49i° 0 ;IJ 

EW n [Fe Ka] (eV) 

163 HI 

E n [Fe K/3] (keV) 

7 034 +0 - 037 
‘ - uo ^-0.038 

In [Fe K/3] (10“° photons cm -2 s _1 ) 

n 43+ 0 - 20 

u.^tO-o 23 

EW n [Fe K/3] (eV) 

97+12 
" ‘ —14 

In [Ni Ka] (10 -5 photons cm' 2 s -1 ) 

< 0.29 

EW n [Ni Ka] (eV) 

< 28 

Reflection fraction, b R(9 obs = 43°) 

<0.60 

Reflection fraction, b R(9 obs ~ 60°) 

<0.78 

*0.5—2 keV ( 10 ~ n er 6 cm_2 s_1 ) 6 

0.16 

*2-10 keV ( 10_11 er g cm “ 2 S- 1 ) 6 

1.15 

£ 0 . 5-2 kev (W 42 erg s' 1 ) 0 

0.70 

^ 2-10 kev (!0 42 erg s _1 ) c 

1.60 


All values and energy ranges in this table (except Fo. 5-2 kev and F 2-10 kev) refer to the rest-frame of NGC 2992 (z = 0.00771). 
Statistical errors are for 90%, one interesting parameter (Ax 2 = 2.706) and were derived with fifteen parameters free, with I n; Fe Ka 
fixed at 0. Upper limits on the intensity and equivalent width (EW) of the Ni K a line were derived with sixteen parameters free. 
° The upper limits on the widths of the line core of Fe Ka and Fe K/3 have been corrected (as described in section 3.2) for the 
degradation of the XIS spectral resolution. b Upper limits on R are for 9 obs fixed at the best-fitting value of 43° and for 0obs = 60°, 
as indicated. c Observed-frame fluxes, corrected for the CCD contamination but not corrected for Galactic and intrinsic absorption. 
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d Intrinsic, rest-frame luminosities, corrected for all absorption components. 



